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Abstract. 

Observations show the increase of high-frequency wave 
power near magnetic network cores and active regions in 
the solar lower atmosphere. This phenomenon can be ex- 
plained by the interaction of acoustic waves with a mag- 
netic field. We consider small-scale, bipolar, magnetic field 
canopy structure near the network cores and active regions 
overlying field-free cylindrical cavities of the photosphere. 
Solving the plasma equations we get the analytical disper- 
sion relation of acoustic oscillations in the field-free cavity 
area. We found that the m = 1 mode, where m is azimuthal 
wave number, cannot be trapped under the canopy due to en- 
ergy leakage upwards. However, higher (m > 2) harmonics 
can be easily trapped leading to the observed acoustic power 
halos under the canopy. 



also been shown that the quiet-Sun chromospheric mag 
netic network elements are surrounded by "magnetic shad 
ows", which lack the oscillatory powe r at highe r fre 



quency range jMclntosh & Judge, 2001; Kr ijger et all 12001 



Vecchio et al 



l2007h . Therefore, both the photospheric 
power halos and the chromospheric magnetic shadows prob- 
ably reflect the same physical process of acoustic wave in- 
teraction with overlying magnetic field (while in subsurface 
regions the rotational and the meridional non-uniform flows 
are supposed to have an impact on the formation of the a cous- 
tic wave power spectra. Shergelashvili & PoedtsL 2005 ). 

The properties of propagating acoustic waves are closely 
related to the magnetic field structure. The numerical 
calculations show that the propagation of acoustic dis- 
turbances in the solar atmosphere is strongly dete r rnined 



by the overlying m agnetic canopy ([Rosenthal et all 12002 : 



Keywords. Sun: photosphere - Sun: acoustic osciUations iBogdan et aU l2003b . The canopy has b een usually mod- 



eled with purely horizontal magnetic field ( Evans & RobertsL 
[1990), but recent high-resolution observatio ns reveal more 



1 Introduction 

Waves play an important role in the dynamics of the so- 
lar atmosphere. Observations show an increase of high- 
frequency power {ly > 5 mHz) in the surroundings of active 
regions in velocity p ower maps sometimes called as photo- 
spheric power halos ( Braun et all 1 992t Brown et al. . 1992; 
Hindman & Brownl 119981 Ijain & Habeii l2002h . The halos 



comp l ex small-scale structu re of the field (iDe Wijn et al. 



boost ICenteno et all l2007h . It has been suggested that 



the magnetic field has small-sc ale closed loop st r ucture 



in the vicinity of netwo rk cores dMcIntosh & Judgel 12001 



Schrijver & TitleL l2003h . The inclined magnetic field may 



were not found in Doppler power maps at lower frequencies 
(3 mHz). Observations also sho w a lack of power halos in 



con tinuum intensity power maps ( Hindm an & Brownl Il998t 
IjaTn & Ha ber. 2002; Muglach et al., 2005). 

On the other hand, recent observations reveal the de- 
crease of the acoustic high frequency power in the chro- 
mosphere and_Jts i ncreas e in the photos p here n ear active 



regions (IMuglachl 120031; iMuglach et all l2005h . It has 
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channel low-frequ ency photospheric osc illations in the chro- 
mosphere/corona dDe Pontieu et al. [ I2OO4I) . 

Here we use a model of small-scale bipolar magnetic 
canopy near a chromospheric network core and/or an active 
region. We suggest that granular cells may form field-free 
cylindrical cavities under the magnetic canopy due to the 
transport of magnetic flux towards boundaries. These cav- 
ities may trap high-frequency acoustic oscillations, while the 
lower-frequency harmonics may propagate upwards in form 
of magneto-acoustic waves. 

Sect. 2 gives the analytical approach, obtained dispersion 
relation and resulting oscillation spectrum. Sect. 3 includes 
discussion and comparison of theoretical findings to obser- 
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1. Schematic picture of the magnetic canopy overlying a field- 



Fig. 

free cavity. 



vations. Sect. 4 briefly summarizes the results. 



2 The model 

We use the ideal magnetohydrodynamic (MHD) equations 
which can be written in the following form 



-=Vx(vxB), 
P-^ +P(v-V)v = 
| + V.(pv).0, 



47r 



(VxB)xB, 



dp 
di 



(v-V)p + 7pV-v = 0, 



(1) 
(2) 
(3) 
(4) 



where v denotes the fluid velocity, B the magnetic field, p 
the pressure, p the mass density, and 7 the ratio of specific 
heats. Gravity effect is omitted from the consideration for the 
sake of simplicity. With this simple model we intend to esti- 
mate the potential of the proposed mechanism. Gravity will 
be taken into account in future developments of the present 
model. 

We assume that the magnetic field is vertical in the 
core of the chromospheric network, but it gets the shape 
of small-scale closed lo o p systems in the surr ounding . 
dMcIntosh & Judgel I2OO1I; ISchriiver & Titlel l2003b . Gran- 
ular motions then may form field-free cavities under the 
canopy as they carry the magnetic flux towards the cell 



boundaries. These field-free cavities have granular dimen- 
sions and may become resonators for acoustic oscillations. 
Consequently, we have two different regions (Fig. 1): 

region I — field-free cavity, r < rg , 
region II — magnetic canopy, r > tq . 

We use a cylindrical coordinate system (r, 0, z) and con- 
sider an unperturbed cylindrical magnetic field Bq in the 
canopy area. The magnetic field has only a component 
which depends on the distance r, i.e. Bq = (0, ^^(r), 0). 
The equilibrium in the canopy is then satisfied if 



= 0, 




where po denotes the unperturbed pressure. To avoid further 
mathematical complications, we consider the unperturbed 
hydrodynamic pressure to be homogeneous. The equilib- 
rium magnetic fie ld is then current-free expressed as follows 
dPiazetaLl [20061) 



^4,0 — ■ 



Equations (l)-(4) are linearized, which yields 
dh 



dt 



= Vx(uxBo), 



(5) 



(6) 



i^(BoV)b+ 



1 

47r 



dpi 
dt 

Pi = 



(bV)Bo, 



PoV-u = 0, 



(7) 

(8) 
(9) 



where b, u, pi and pi are the perturbations of magnetic field, 
velocity, pressure and mass density respectively, while pQ is 
a uniform unperturbed mass density and cq = (7Po/po)^^^ 
corresponds to the homogeneous sound speed. Note, that a 
background flow is absent in our consideration. 

For simplicity, we consider the 2D case and restrict the 
analysis to the (r, 0)— plane. In principle, the z— direction 
can also be considered, but this further complicates the pre- 
sentation (namely, resonant absorption may take place) and it 
is left for future considerations. Then r— and (/(—components 
of equations (6)-(9) are given by 



dbr Bcj, dur 
dt r dd) ' 



db^ 
dt 



-B, 



Po- 



dUr 

dT 



dUr 

' dr 
2dpi 



B, 



Ur 
> ; 

r 

B^ db^ 
4tt dr 



B^ dbr 
Aur dcf) 



(10) 



(11) 



(12) 
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dt 



Cq dpi 
r dd> 



dpi 
dt 



Pn- 



dUr 

dr 



Ur 

Po — 



Po du^ 



= 0. 



(13) 



(14) 



To get the oscillation spectrum in the cavity area, we have 
to solve the equations in the cavity and canopy regions sepa- 
rately and then merge the obtained solutions at the interface 

(r = ro). 

Medium is field-free in the cavity (r < ro) and therefore 
can be described by pure hydrodynamics. Then equations 
(12)-(14) lead to the Bessel equation after Fourier analysis 
with respect to both the t (time) and (p coordinates (the mag- 
netic field is set to zero): 



d'pi 

Qj.2 



idpi 

r dr 



Pi =0, 



(15) 



where uj is the wave frequency, and m is the azimuthal wave 
number 

In the magnetized canopy region (r > ro), equations (10)- 
(14) lead to the Hain-Lust equation 

d_ 

dr 



2 9 



1^2 _ 



dr 



d_ 

dr 



9 2 

m^Cn 



/r2 



Ur 

r 



Ur = 0, 



(16) 



where va = B^/^/A¥po denotes the Alfven speed. 

Analytical solution of equation (16) is complicated. 
Therefore, for simplicity, we suppose that the magnetic en- 
ergy is much higher than the hydrodynamic one within the 
canopy area. Consequently, we use the zero /3 approxima- 
tion hereinafter. The unperturbed configuration must be in 
equilibrium, therefore the unperturbed hydrodynamic pres- 
sure of the cavity and the magnetic pressure of the canopy 
must be balanced at the interface (r = ro) 



CpiPOl 
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where 

coi = 



Stt ' 



Poi '' 



(17) 



(18) 



and coi and poi are the sound speed and the plasma mass 
density in the cavity. Equation (17) then gives the relation 
between the sound and Alfven speeds at the interface 



/7 P02 
Coi = W X VA2, 

V 2 poi 

where 



VA2 



B^{ro) 



(19) 



(20) 



is the Alfven speed at the interface and po2 is the plasma 
mass density in the canopy. 



2.1 Analytical Solutions 

In the field-free cavity area under the canopy (region I), there 
are only acoustic waves. The r-dependence of the mass den- 
sity perturbations in the acoustic waves is governed by equa- 
tion (15) which has the general solution 



Pi = ciJmikir) + C2Ym{kir), 



(21) 



where Jm{kir) iYm{kir)) is the bessel function of the first 
(second) kind. 



fci = — 

coi 



(22) 



and ci, C2 are arbitrary constants. The perturbation must be 
finite at r = 0. Therefore, C2 — and only the first term on 
the right hand side of expression (21) is non- vanishing. 

The radial velocity component of the acoustic wave in the 
region I can be obtained using equations (12) and (21) as 
(setting the magnetic field to be zero) 



Url 



-cij;„(fcir) 



(23) 



where the prime denotes the derivation with respect to r. 

In the magnetic canopy (region II), where the cold plasma 
approximation is used, equation (16) gives 



92 



Qj,2 



1 dilr 

r dr 



1 



v\{r) 



= 0. 



This equation can be rewritten as 



ld_ 

r dr 

where 



d^p 
dr 



m 



21 



Vj^{r) 

= rip, v\{r) 



^ = 0, 



(24) 



(25) 



B 



00 



47rpo2 



The solutions o f this equation are the Bessel f unctions of 
half i nteger order ( Abramowitz & StegunL 1967 : Dfaz et al 



20061) . In region II, only the outgoing wave is physically 
appropriate. Therefore, we choose the Hankel function 

Ur2 = iC'irHra/2{k2r), (26) 

where C3 is an arbitrary constant and 



k. = 



2WA2 



(27) 



The total pressure perturbations in the regions I and II are, 
respectively. 



Pi = CgiPl, 

and 

B^b^, _ ipo2VA2 



4tt 



dii. 



t2 



dr 



Ur2_ 

r 



(28) 



(29) 



Thus the expressions (23), (26), (28), and (29) give the 
transverse velocity and the total pressure perturbations in the 
considered two regions. 
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Fig. 2. Temporal dynamics of velocity perturbations in the field-free 
cavity areas for m — 1,2,3 harmonics. 



Fig. 3. The ratio of imaginary uji and real ujr parts of wave fre- 
quency vs the ratio of Alfven and sound speeds for m = 1, 2, 3 
harmonics. 



2.2 Dispersion relation 

The continuity of the velocity and the total pressure pertur- 
bations at the interface (r — tq) leads to 



CoiPi 



du 



r2 



dr 



Ur2 

r 



Url — Ur2- 



(30) 



(31) 



The substitution of the expressions pi, Uri,Ur2 into equa- 
tions (30) and (31) then gives 

cliCiJjnikiro) = —C3rok2H'^,r^{k2ro), (32) 



cikiJ^{kiro) ^ czroH^/2{k2ro). (33) 



The condition for a non-trivial solution of equations (32)- 
(33), then yields the general dispersion relation, viz. 



via P02 kiJI^{kira) 



fc2g^,/2(fc2ro) 
Hm/2{k2ro) 



(34) 



The dispersion relation (34) is a transcendental equation 
for the complex w. The imaginary part of wave frequency lu 
indicates a wave leakage from the field-free cavity into the 
ambient magnetic canopy. The question then arises how im- 
portant this wave leakage is and which are the most sensitive 
parameters for trapping (or leakage) of the waves under the 
canopy. The analytical solution of equation (34) is compli- 
cated. Therefore, we apply numerical techniques to solve it. 

The properties of the dispersion relation (34) depend on 
the azimuthal wave number m and the ratio between sound 



(coi) and Alfven (va2) speeds (or the ratio between poi and 
P02, see Eq. (19)). The dimension of the field-free cavity, 
ro, also stands as a free parameter, but it influences only the 
wave periods. 

Numerical solution of equation (34) shows that the ratio of 
imaginary uji and real uJr parts significantly depends on the 
azimuthal wave number m. The ratio is higher for m — 1 and 
quickly decreases with increasing m. Fig. 2 shows the tem- 
poral dynamics of radial velocity perturbations in the field- 
free cavity region for modes corresponding to different m- 
values. The amplitude of the first (m = 1) harmonic quickly 
decreases, which indicates the rapid radiation of this mode 
into the overlying canopy. On the other hand, the m — 2 
mode undergoes a very small leakage and the m — 3 mode 
has almost no leakage. Thus, the first harmonic of the acous- 
tic oscillations cannot be trapped in the cavity area due to the 
rapid leakage into the canopy, while the higher harmonics, 
i.e. those with to > 2, can be easily trapped, which may lead 
to the observed increased acoustic power. 

Fig. 3 shows the ratio of the imaginary uji and real ojr parts 
vs the ratio of the Alfven and sound speeds for different har- 
monics. We see that the decrease of the ratio between the 
Alfven and sound speeds leads to an enhanced leakage. 

Fig. 4 shows the dependence of real (top panel) and imag- 
inary (bottom panel) parts of the frequency on the size tq of 
the field-free cavity region under the canopy. It is evident 
that the real part of the frequency decreases with increasing 
ro, as can be expected from physical considerations. 
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3 Discussion 



We have studied the spectrum of acoustic oscillations in a 
cylindrical field-free cavity under a small-scale bipolar mag- 
netic canopy in the solar atmosphere. It is shown that the 
m = 1 (where rn is the azimuthal wave number) harmonic 
of the acoustic oscillations cannot be trapped in the cav- 
ity as a result of the energy leakage in the upward direc- 
tion. The energy radiation occurs through the propagation of 
fast magneto-acoustic waves in overlying magnetic canopy. 
However, higher m> 2 harmonics can be trapped in the cav- 
ity, leading to the observed increased high-frequency power 
in the photosphere. 

There are three different explanations of power halos 
proposed in the literature: (i) the enhance ment of acous 



tic emission by some unknown source (Braun et al.L 11992 



Brown et al.lll992l:|jain & Habeii l2002h . (ii) incompressible 



oscillations, such as Alfven waves or transverse kink waves, 
in magnetic tubes (Hindman & Brown, 1998.) . and (iii) the 
interacti on of acoustic waves with the overlying magnetic 
canopy ("Muglac h et al.[l2005l) . 

We suggest that the surroundings of magnetic net- 
work cores and active regions consi st of many small- 
scale closed magnetic canopy structures ("Mcintosh & Judge}, 
I2OO I: Schrijver & Title, 2003). Granular motions trans- 
port the magnetic field at boundaries and consequently cre- 
ate field-free cylindrical cavity areas under the canopy (see 
Fig. 1). These field-free cavities may be filled by trapped 
acous tic oscillations stochastically excit ed by granular mo- 
tions (lLighthill[[l952llSelwa et al.Ll2004 . For a typical pho- 
tospheric sound speed of ~ 7 km/s and a typical granular 
radius of ~ 400 km, the period of the m — 2 harmonic is 
~ 3 min and the period of the m — 3 harmonic is ^ 2 min. 
The range of observed enhanced power frequency is about 
5 — 7 mHz (with periods of 2 — 3 min), which is in agree- 
ment with our findings. The m = 1 harmonic is leaky and, 
therefore, can not be trapped under the canopy. 

If power halos are related to the acoustic waves then ob- 
servations should show an enhancement in continuum inten- 
sity power maps as well, but in op posite observations s how a 
lack of power halos i n the maps ( Hindman & Brownl Il998 



Jain & HabeA 120021; iMuglach et al.L 12005 '). what needs an 
adequate explanation. Here, we suggest a natural explanation 
for this phenomenon. Indeed, the antinodes of mass density 
and velocity components in standing acoustic oscillations are 
located at different places: the location of the maximal ve- 
locity oscillations corresponds to the location of the zero 
amplitude intensity oscillations and vice versa. Therefore, 
the enhanced Doppler velocity power at fixed height of the 
atmosphere automatically suggest the absence of intensity 
power at the same height. This suggestion can be checked 
by searching intensity power halos at different heights from 
the surface. We believe that new observations with high res- 
olution from the Hinode spacecraft will shed light on this 
problem. 
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Fig. 4. Real Ur (upper panel) and imaginary uji (lower panel) parts 
of wave frequency vs the radius of field-free cavity ro for m = 



However, the formation of power halos due to incom- 
pressible oscillations can not be completely ruled out. 
Recent theoretical and two-dimensional numerical sim- 
ulations outlined the importance of the /? ~ 1 region in 
the solar atmosphere, where wave conversion or reflec 
tion occurs feosentha l'erall l2002t bogdanetall |2003| 



Zag arashvih & Robertsl 2006; Kuridze & ZaqarashviU , 



i2007l) . Indeed, our consideration implies that the hydrody- 
namic and magnetic pressures have approximately the same 
value at the interface (see Eq. 17). Therefore, the acoustic 
waves may transfer energy into th e incompressible waves 
through non-linear interactions |Zaaarashvih & Robertsl 



20061 Kuridze & ZaqarashviU . 20071) . which may lead to the 
observed power halos. However, the observed oscillations 
hardly show non-linear behavior, which complicates the 
explanation of power halos by incompressible oscillations 
considerably. 

Another important alternative for the proposed mecha- 
nism can be r esonant ab s orption of waves in inhomoge- 
neo us plasma dTirrv et all 1 19981; IPinter & Goossensl Il999t 



Pint er et alil2007h . This process can also be particularly im- 



portant for 3D consideration of proposed model, when one 
considers the wave propagation in the z direction as well. 
This is out of the scope of present paper, but would be inter- 
esting to study in the future. 

It must be mentioned that the equilibrium used in this pa- 
per is simplified as gravitational stratification, which is im- 
portant in the solar atmosphere, is ignored. The stratifica- 
tion leads to a Klein-Gordon equation for the propagatin 
waves with a cut-off for wave frequencies dRoberts. ,200- 



ng 
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Erdelyi et al.Ll2007l) . However, wave propagation in inclined 
magnetic field may lead to decrease of cut-off frequencies 
and simplifies the penetration of low er- frequency oscillations 
in higher regions (IDe Pontieu et al.l 12004). Non of the har- 
monics has purely vertical propagation in our model because 
of the non-zero azimuthal wave number m; all harmonics 
propagate with an angle to the vertical (note, that all har- 
monics have standing wave behavior along azimuthal direc- 
tion). Therefore their frequencies are above the cut-off value 
and consequently the harmonics are not evanescent. How- 
ever, inclusion of gravitational stratification is necessary for 
a more profound understanding of the wave trapping in cav- 
ities. This wiU be the next step of our study. 



4 Conclusion 

We have studied the spectrum of acoustic oscillations in 
the cylindrical field-free cavity regions under the small-scale 
magnetic canopy near the magnetic network cores and active 
regions. We found that the first harmonic of acoustic oscilla- 
tions cannot be trapped in the cavity due to the energy radia- 
tion by fast magneto-acoustic waves in the canopy. However, 
the higher (m > 2) harmonics can be trapped there, lead- 
ing to the observed enhancement of high-frequency acoustic 
power in the photosphere. Future detailed study of the pro- 
posed mechanism including gravitational stratification and 
3D models is necessary. 
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